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1 Introduction 

In this short communication, I give a very simple derivation of the Jarzynski 
equahty P 

Ai^ = -r'ln(e-^^), (1) 

which aUows to compute the free energy difference AF = Fb — Fa of a body, 
which is driven between two equihbrium states A and B by an external (time- 
dependent) force A, from the probability distribution function P(M^) of the work 
done on the system 

regardless of the nature of the transformation (reversible or irreversible) between 
the states A and B. Here Ht{x,X) denotes the full (time-dependent) Hamilto- 
nian of the system, which is supposed to be in thermal contact with a heat bath 
of temperature T = l/ksP, the whole ensemble being isolated, x = {p,q) = 
(pi, . . . ,Pn, , ■ ■ ■ 1 q") denotes the phase-space point of the system (which is 
supposed to have 2n degrees of freedom), and {f{W)) = J f{W) P{W) dW de- 
notes an average over the ensemble of realizations of the work W done on the 
system (for any "smooth" function /). 

There are two parts in this derivation. The first is mathematical and will be 
discussed later. The result is the following: let £'(2;) be the (time-independent) 
energy of a body in thermal equilibrium with a heat bath of temperature T. 
If its energy can be written in the additive manner E{x) ~ Eq{x) + V{x), 
which can always be achieved by making certain very general assumptions, then 
one can compute the free energy F of the body (and all the relevant ther- 
modynamic quantities) from the only knowledge of the equilibrium Gibbsian 
distribution ^0(2;) = Zq^ e~^^°^^\ where Zq = J e~^^°^^^ dx is the partition 
function corresponding to the "unperturbed" energy Eq (x) of the body (and to 
its "unperturbed" free energy _Fo = ^P^^ biZo). 
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The second part of the argument, which relates the perturbation energy V to 
the thermodynamical work W done on the system, first in the time-indcpcndcnt 
case, and then in the time-dependent case, is physical. 

2 Thermodynamic perturbation theory 

The fundamental need for perturbation theory runs as follows: for example, 
in quantum mechanics, the goal is to compute the eigenvalues A (the energy 
spectrum) and the eigenvectors u (the wave functions) of the (time-dependent 
or not) Hamiltonian H{x) of a system, from the eigenvalue equation H{u) = 
A u, where H is generally a very complicated differential operator, occasionally 
depending on time via some time-dependent energy coupling (due, e.g., to some 
external time-dependent force). 

However, for evident analytical reasons, the calculus of the eigenvalues A and 
of the eigenvectors u is generally intractable. In order to circumvent this prob- 
lem, one generally write (when this is possible) the energy of the system in the 
additive manner H{x) = Ho{x) + V{x), where Ho{x) denotes the "unperturbed" 
energy of the system (i.e. the quantity from which one can compute the eigen- 
values Ao and the eigenvectors uo), and V{x) the perturbation energy. Then, 
one computes by iteration A and u from the only knowledge of the "unperturbed" 
eigenvalues Aq and eigenfunctions uq, which are computable Eventually, the 
perturbation energy V{x) can depend on some external force A, time-dependent 
or not (e.g., A can be an electrical or magnetic field). 

The situation is very similar in statistical mechanics: once ones has set 
up the Hamiltonian E{x) of the system (which, again, is in thermal contact 
with a heat bath of temperature T, the whole ensemble being isolated), the 
fundamental problem is to compute the partition funtion Z of the system 




where dx = dpdq ^ Yi ^Pi , from which one can compute the free energy F 
of the system 

F -f3-HnZ, 

and, e.g., the thermal average energy of the system (E) = — ^^^^ 

Again, these well-known formulas are purely formal, and inapplicable in real- 
istic cases, because the energy E(x) of the body is generally a very complicated 
function of the coordinates q and momenta p of the "particles" constituting the 
body, and the integration over the phase-space cannot be carried out in order 
to compute Z [in the quantum case, this calculus is of combinatoric nature, 
and requires a counting of the (occasionally degenerate) states |n) of the system 
via the formula Z = X]^"'''^", which also becomes quickly an impossible task]. 
However, as we will show below, much more can be said about Z, F and (E), 
at least formally. 
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Let us first assume that the energy of the body (which, again, is in thermal 
contact with a heat bath of temperature T, the whole ensemble being isolated) 
is time-independent, and that it can be written in the additive manner 

E{x)^Eo{x) + V{x), (3) 

where "everything" can be calculated from the "unperturbed" energy Eq{x) of 
the body [the partition function Zq = J e"^^"^^'' dx, the free energy Fa = 
— /3^^1nZo, etc.], but not from its "full" energy E{x). Therefore, the funda- 
mental problem is to devise a method to compute, even approximatively, the 
partition function Z, the free energy F, etc., corresponding to the "full" energy 
E{x) of the system. Indeed, wc will show that if Eq. (j2Jl holds, then one can 
compute exactly Z and F, without anymore assumptions than Eq. (j2J), which 
relies on the separability of the Hamiltonian of the system, and on thermal 
equilibrium. 

From these very few assumptions, the algebra is very simple. Let us first 
compute the partition function Z = J Q-P^i^) da; of the system. From Eq. Q, 
one gets 

Z = Zaj po{x) e-''^^-) dx = Zo (c-''^)„ , (4) 

where po{x) = Zq"^ c"^^"^^' is the "unperturbed" equilibrium Gibbsian distri- 
bution of the system. Therefore, the free energy of the system is given by 

F = ~/3-HnZ = F„-fr^ ln(c-^'')o- (5) 

This result was derived up to the second order [i.e. only with the first two 
cumulants of V in Eq. ((SJ] by Peierls and Bogoliubov, who gave an upper bound 
for the free energy F of the system P] [indeed, Eq. ^ indicates that the "full" 
free energy F is merely the cumulants generating function of the perturbation 
energy V]. Our result is much stronger, and looks amazingly similar to the 
Jarzynski equality JQ). 

3 Dynamical interpretation 

The preceeding result holds for a time-independent and equilibrium system. 
However, provided that the system is initially (say, at time t = 0) at equilibrium, 
these results are also true if the energy of the system depend on time through 
a time-dependent perturbation energy, which can be due to an external force 
A which is switched on at times t > 0. Therefore, if the energy of the system 
reads 

Et{x,X)^Eoix) + Vt{x,X), (6) 

where the time-dependent perturbation energy Vt is only non-zero for times 
t > 0, then one trivially has 

Ft = Fo~ In (c-'^^*)^ for all t > 0, (7) 
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which has exactly the same meaning as Eqs. (@J-(jSl- 

In fact, the perturbation energy Vj can be identified to the work Wt done 
on the system on the time interval (0, t) in a very simple way [one simply has 
to put tA = and ts = t in Eq. (|2J), but here B does not necessarily denotes 
an equilibrium state]. Since the energy Et of the system depends explicitly on 
time only through the external force X{t), one has by applying the chain rule 

dEt{x,X) _ dEt{x,X) _ dE{x, A) dA 
dt ^ dt ^ dX dt ■ 

Therefore, by Eq. Q one has Wt ^ Et — Eq ^ Vt, and one finally gets 

^^, = Fo-rMn(e-^^*)g, (8) 
which resembles amazingly the Jarzynski equality l|T]l. 

4 Discussion and conclusion 

The formula © obtained above is not really the Jarzynski equality ^ , since the 
Jarzynski equality computes the free energy difference of the system from the 
work distribution P(M^), which takes into account all the (in principle) noncqui- 
librium fluctuations of the work W done on the system during the transfor- 
mation between the (equilibrium) states A and B. In our derivation, this is 
obviously not the case, since the average (■) is performed from the equilibrium 
Gibbsian distribution pq{x) = Z^^ c-l^^°^''\ 

However, this slight detail perhaps throws some light on the debate started 
by Cohen and coworkers ^ , and perhaps explains why all the experiments which 
has been performed in order to test the Jarzynski equality are successful, even 
in very defavorable (irreversible) cases (and more strikingly, when the state B 
is not an equilibrium one) 

To conclude, we should add that among the hypotheses needed to the result 
the crgodicity of the system is of course required. However, even unrealistic, 
this is a rather weak statement, because we only need crgodicity at times t < 0, 
that is to say without external (time-dependent) driving force A. 
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